THREE  DIMENSIONAL  STRESS  ANALYSIS 
OF  AN  ARTILLERY  PROJECTILE  JOINT 
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ABSTRACT 

“^This  study  presents  the  results  of  a three-dimensional  elastic-  p 
plastic  dynamic  stress  analysis  of  one  of  the  structural  joints  encoun- 
tered in  artillery  projectiles.  The  particular  spline  joint  analyzed 
has  equally  spaced  set  screws  around  the  surface  of  th-  projectile.  In 
general,  these  types  of  structural  joint  problems  are  variable  contact 
problems  in  that  the  interaction  between  the  set  screws  and  their  bearing 
surface  along  with  the  interaction  between  the  interfaces  of  the  joint 
are  nonlinear  in  nature.  Due  to  the  complexity  of  the  structural 
configuration  and  loadings  of  the  joint,  the  finite  element  method  has 
been  used  to  solve  the  problem.  The  numerical  analysis  covers  the  time 
from  initial  launching  to  barrel  exit.  Stresses  and  deformations  m tne 
joint  are  determined  at  various  stages  of  loading.  The  effect  of  the  set 
screws  and  set  screw  holes  on  the  stress  distributions  in  the  joint  is 
examined  in  detail.  ^ 

INTRODUCTION: 

Artillery  projectiles  are  subjected  to  extremely  high  loads  during 
firing.  At  present,  the  design  of  artillery  projectiles  is  greatly 
facilitated  by  the  exploitation  of  the  finite  element  method.  But  the 
application  of  the  method  has  been  limited  to  the  simplified  two  dimen- 
sional or  axisymmetric  analysis.  This  is  due  to  the  complexity  of  the 
geometries  of  the  projectiles  thus  requiring  very  long  computer  analysis. 
Although  successes  were  achieved  in  the  previous  designs,  there  have 
always  been  concerns  of  the  structural  integrity  of  the  projectiles  viith 
the  presence  of  high  local  stresses.  In  the  case  of  XM753  projectiles, 
high  local  stresses  appear  in  the  region  of  the  pinned  joint.  In  the 
case  of  XM783  projectiles  the  use  of  spline  joint  with  set  screws  in  the 
design  also  leads  to  high  local  stresses  in  the  joint  area. 

After  firing,  an  artillery  projectile  is  subjected  to  various  loads 
continuously  changing  with  time.  It  experiences  first  a very  high 
compressive  load  in  the  axial  direction  during  in  bore  flight  and  then  a 
high  tensile  load  at  barrel  exit.  Such  a loading  history  causes  the 
projectile  to  undergo  a stress  reversal  i.e.  from  a stress  state  of 
compression  to  that  of  tension.  However  it  has  been  a common  practice  in 
the  design  of  artillery  projectiles  to  perform  either  one  or  two  simple 
independent  stress  analysis.  The  one  which  is  generally  carried  out  is 
a two  dimensional  quasi-static  analysis  in  which  the  loads  at  the  peak 
linear  acceleration  during  in  bore  flight  are  used  to  compute  the 
stresses  and  deformations  in  the  projectile. 


Sometimes  a two  dimensional  dynamic  analysis  is  also  conducted  to 
determine  the  stresses  and  deformations  in  the  projectile  caused  by 
the  sudden  drop  off  of  pressure  of  the  propellant  gas  at  barrel  exit. 

In  the  present  investigation  a detailed  three  dimensional  elastic- 
plastic  dynamic  stress  analysis  of  a spline  joint  of  XM785  projectile 
is  performed.  The  purposes  of  the  study  are.- 

(1)  to  verify  the  structural  adequacy  of  a proposed  design. 

(2)  to  asses  the  effect  of  set  screws  and  set  screw  holes  on  the 
stress  distributions  in  the  joint  region  of  the  projectile. 

In  addition  the  present  work  also  serves  as  an  initial  effort  to 
determine  the  extent  of  the  influence  of  the  stress  reversal  or 
Bauschinger  effect  on  the  structural  integrity  of  the  projectile. 

DESCRIPTION  OF  LOADING  CONDITIONS 

After  a projectile  is  fired  and  before  it  departs  the  gun  barrel 
it  is  subjected  to  a combination  of  the  following  loads: 

(1)  axial  compressive  load  due  to  linear  acceleration  of  the 
pro j ect i 1 e 

(2)  centrifugal  load  due  to  angular  rotation  of  the  projectile. 

(3)  torsional  load  due  to  angular  acceleration  of  the  projectile. 

(4)  internal  load  due  to  interaction  of  interior  components  and 
projectile. 

(5)  external  load  due  to  gun  tube  constraint,  rotating  band  pressure 
and  balloting. 

As  the  projectile  departs  the  gun  barrel,  it  experiences  a tensile  load 
or  a negative  set-back  load  (elastic  release)  in  the  axial  direction 
resulting  from  the  sudden  drop  off  of  propellant  pressure  at  the  barrel 
exit.  Among  the  loads  the  axial  load  is  the  dominating  one. 

It  has  been  a general  practice  to  omit  the  effect  of  the  torsional 
load  induced  by  the  angular  acceleration  in  the  design  of  the  projectiles. 
Analysis  has  shown  that  as  a result  of  such  an  omission,  the  magnitude  of 
the  effective  stresses  in  the  projectiles  are  about  2%  - 4%  lower  (1). 
Since  the  present  investigation  concerns  the  determination  of  the  stress 
distributions  in  the  region  of  the  joint,  only  the  portion  of  the 
projectile  in  the  neighborhood  of  the  joint  is  considered.  The  area  where 
the  rotating  band  is  located  is  excluded.  Thus  the  load  due  to  rotating 
band  pressure  is  not  included  in  the  analysis.  Only  the  axial,  centrifugal 
and  internal  loads  during  in  bore  flight  and  the  negative  set  back  load  at 
barrel  exit  are  considered  in  the  three  dimensional  dynamic  stress 
analysis  of  the  joint.  Fig.  1 shows  the  linear  acceleration  of  the 
projectile  used  in  this  analysis  for  the  calculation  of  the  axial  load. 

The  projectile  reaches  an  acceleration  of  1 7,000g  in  about  6 mi  1 1 i seconds , 
zero  acceleration  at  barrel  exit  and  immediately  is  subjected  to  a 
negative  acceleration  of  2000g. 
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METHOD  OF  STRESS  ANALYSIS 


Figs.  2 and  3 shows  respectively  the  geometric  configuration  of 
an  artillery  projectile  and  a typical  section  of  its  joint.  The  latter 
is  formed  by  passing  a plane  at  the  midpoint  of  two  neighboring  set 
screws  and  another  plane  through  the  center  of  one  of  the  set  screws. 

The  two  planes  are  parallel  to  the  axis  of  the  projectile. 

In  the  initial  phase  of  the  analysis,  the  finite  element  model  of 
the  joint  (see  Fig.  4)  was  created  by  using  the  computer  program  PATRAN-G 
(2).  It  employs  sophisticated  interactive  color  graphics  and  a powerful 
geometry-based  language  for  geometry  construction  and  finite  element 
modeling  during  the  first  phase  of  the  program.  Then  nodal  point  and 
finite  element  generation,  assignment  of  physical  properties,  application 
of  external  loads  and  definition  of  sliding  interface^  are  accomplished 
in  the  second  phase  of  the  program. 
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The  explicit  three  dimensional  finite  element  code  DYNA3D  (5)  was 
used  to  compute  the  stresses  and  deformations  in  the  finite  element 
model  of  the  joint.  DYNA3D  is  designed  to  analyze  the  large  deformation 
dynamic  response  of  inelastic  solids.  It  has  a contact  algorithm  that 
can  model  gaps  and  sliding  materials  interfaces.  It  uses  a 8-node  constant 
stress  solid  element  and  one  point  integration  in-  element  stiffness 
calculations.  It  is  programmed  to  take  full  advantage  of  vector 
optimization  on  the  CRAY-1  (a  class  VI  machine)  and  can  execute  at  less 
than  0.67CPU  (central  processor  Units)  minutes  per  million  mesh  cycles. 

A symmetric,  penalty  based,  contact- impact  algorithm  was  implemented 
that  not  only  reduced  hourglassing  problems,  but  also  was  considerably 
faster  in  execution  speed  and  exceedingly  reliable. 
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There  were  1773  eight  node  brick  elements  and  2611  nodes  in  the 
finite  element  model  of  the  joint.  Four  sets  of  gap  or  sliding  interfaces 
were  required  in  the  model.  Three  sets  of  sliding  interfaces  were 
used  to  model  the  contact  between  the  case  structure  and  the  wedge.  A 
sliding  interface  with  small  initial  gap  was  used  to  model  the  interaction 
between  the  set  screw  and  the  hole. 

The  computations  were  performed  on  a CRAY-1  computer  at  Sandia 
National  Laboratories.  The  computer  time  required  to  complete  the  dynamic 
analysis  was  about  2.8  CPU  hours. 

DISCUSSION  OF  RESULTS 
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The  post  processor  program  GRAPE  (4)  was  used  to  obtain  plots  of 
stress  contours  on  the  surfaces  of  the  joint  and  the  deformed  shapes  of  the 
joint  at  the  time  when  the  maximum  linear  acceleration  is  reached  and 
tlso  at  barrel  exit.  These  are  shown  in  Figs.  5 to  11.  Examination  of 
the  results  of  the  analysis  leads  to  the  following  findings: 

(1)  Due  to  the  presence  of  set  screws,  a larger  portion  of 

the  axial  compressive  load  is  transmitted,  during  in-bore 
flight,  through  the  center  sections  of  the  regions  bounded 
by  each  pair  of  set  screws.  Consequently,  the  center 


E; 


sections  experience  higher  stresses.  This  is  different 
from  the  uniform  stress  distribution  in  the  circ..  fer ,ntial 
direction  found  in  the  two  dimensional  stress  analysis. 

2.  Higher  stresses  occur  on  the  outer  surfaces  of  the  case 
structure  and  the  wedge  at  both  peak  axial  cojt ression 
and  tension  at  barrel  exit  due  to  bending  effect. 

3.  Local  yielding  occurs  in  the  set  screws,  s^t.  .crew  holes 
and  other  areas  in  the  case  structure  and  t'.e  wedge  at 
time  of  peak  axial  compression  and  that  of  tension  at 
barrel  exit.-  Initial  examination  of  the  pia'tic  strains 
(or  Bauschinger  effect)  at  peak  axial  comp: r ision  leads 
one  to  believe  that  they  have  a negligible  effect  on  the 
level  of  stresses  in  the  joint  at  barrel  f .it.  This  is 
because  their  magnitudes  are  small  (.001)  ..ind  their 
locations  are  also  different  from  the  pla  ;ic  strains 
induced  by  the  tensile  load  or  negative  s’-t-back  at 
barrel  exit. 

4.  Similar  stress  distributions  were  obtai'c.  when  a wedge 
made  of  either  steel  or  titanium  was  used  with  a titanium 
case  structure. 

In  summary,  a detailed  three  dimensional  e\  ist i c-pl ast i c dynamic 
stress  analysis  of  an  artillery  projectile  was  bxained  to  determine 
the  adequacy  of  a proposed  design.  The  analysis,  was  performed  by 
using  the  DYNA3D  finite  element  program.  The  results  of  the  present 
analysis  will  be  verified  by  laboratory  test  currently  under  preparation 
at  AMMRC.  Future  analytical  work  will  include  application  of  fracture 
mechanics  to  the  problem  and  accurate  determination  of  the  Bauschinger 
effect  by  using  actual  stress  strain  curves  of  the  materials.  The 
present  analysis  employed  isotropic  strain  hardening  which  does  not 
accurately  represent  the  material  behavior  in  a loading  condition  where 
there  is  a reversal  of  applied  loads. 
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See  Table  1 for  Values  of  Symbols 
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